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When two cavitation bubbles exist in a confined space, the interaction between the bubbles significantly affects 
the characteristics of bubble dynamic behaviors. In this paper, a three-dimensional (3D) model is established to 
study the growth and collapse of two cavitation bubbles in a heated tube and its effects on heat transfer. The liquid 
and gas phases throughout the calculation domain are solved by a set of Navier-Stokes equations. It is assumed 
that the gas inside the bubble is compressible vapor, and the surrounding liquid is incompressible water. The mass 
transfer between two phases is ignored. The calculated bubble profiles were compared to the available experi- 
mental data, and a good agreement has been achieved. Then, the relationship among the bubble motion, flow field 
and pressure distributions was analyzed. On this basis, the effects of bubble interaction on the heat transfer be- 
tween the wall surface and sounding liquid were discussed. It is found that heat transfer in the centre wall region 
is enhanced owing to the vortex flow and micro-jet induced by the bubble contraction and collapse. In contrast, 
the highest surface temperature appears in the surrounding region, which is mainly attributed to the thermal resis- 


tance induced by the bubble. The present study is helpful to understand the heat transfer phenomenon with cavita- 


tion in the liquid. 
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Introduction 


The existing of caviation bubbles in a liquid can cause 
noise and erosion, which is very harmful and should be 
avoided in the hydrodynamic system!'~!, However, on the 
other hand, those bubbles can also be widely applied in 
many fields, such as ultrasonic cleaning, waste water 
degradation, shock wave lithotripsy, and enhancing heat 
transfer", Therefore, extensive studies on the cavita- 
tion bubble dynamics have been conducted to reveal the 
mechanism of cavitation. The development of high speed 
photography technique makes it possible to capture the 
extremely complex dynamic processes of cavitation bub- 
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bles. Brujan et al.!'7! investigated the behavior of the la- 
ser-induced cavitation bubbles near a boundary. They 
pointed out that the shape of the liquid jet formed after 
the bubble collapse is strongly dependent on the dimen- 
sionless distance between the bubble and boundary. Then, 
the penetrating jet and its impacts on the opposite wall of 
the non-spherically bubble collapse were further studied. 
The maximum jet velocity and shock pressure up to 280 
m/s and 25 GPa during the collapse stage have been ob- 
served''?!, Yang et al.''*! employed the hydrophone sys- 
tem with a sample rate of 1 billion samples per second to 
monitor the dynamic features of a cavitation bubble near 
the solid boundary. The energy loss of the bubble from 
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Nomenclature 

total energy (J) 

force (N) 

heat transfer coefficient (W/m’-K) 
distance from the centre point to the bubble 
effective conductivity (W/ m-K) 
length of the tube (mm) 
molecular weight (kg/mol) 

static pressure (Pa) 

heat flux (W/cm?) 

bubble radius (mm) 

gas constant (J/K-mol) 

time (s) 

temperature (K) 

velocity (m/s) 

x coordinate (m) 
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y coordinate (m) 


the first maximum expansion to the second maximum 
expansion has been estimated. The cavitation bubble in a 
liquid gap and a confined space with more than one 
boundary has also been widely studied. It is found that 
boundary layers will hinder the bubble collapse as well 
as the formation of the micro-jet!'*!, Due to the short du- 
ration of bubble evolution, most high speed cameras are 
not fast enough to capture the actual time of bubble dy- 
namics. Numerical simulation provides a powerful tool 
that can more accurately reveal the mechanism of cavita- 
tion by setting the parameter ranges that are difficult to 
achieve in the experiment. The Rayleigh-Plesset equation 
was proposed to describe the shape evolution of the bub- 
ble and liquid jets'’®!, Furthermore, the theory was mod- 
ified by considering the contact angle, surface tension 
and kinematic viscosity, which makes the results agree- 
ing better with the experimental observations!'”!, With the 
development of numerical techniques, more methods 
have been widely applied on cavitation bubble dynamics, 
such as the level-set, BEM, Front tracking and VOF!'*?", 
Researches mentioned above lay a solid foundation for 
the application of cavitation. In recent years, it is found 
that the cavitating flow can significantly enhance heat 
transfer, which has drawn a wide attention from scho- 
lars””, Liu et al.'??! investigated the behavior of a single 
bubble near a heated boundary. Results show that the 
micro-jet formed during the bubble collapse is mainly 
responsible for the reinforcement of heat transfer. On this 
basis, the single bubble in a narrow gap and a micro-tube 
with constant heat fluxes are continually investigated to 
reveal the mechanism of the heat transfer enhancement 
by cavitation“. Due to the complexity of the bubble 
dynamic model, most studies are on the simplifying as- 
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Z z coordinate (m) 

Greek letters 

pP Density (kg/m’) 

a volume fraction (-) 

u viscosity (kg/m:s) 

k curvature of the interface (1/m) 
surface tension (N/m) 

Subscripts 

f fluid 

g gas phase 

1 liquid phase 

w wall surface 

max maximum value 

t total 

0 initial 

00 far field 


sumption that the bubble is in a relatively isolated space 
without considering the interaction effect. However, 
when the cavitation occurs in the liquid, the bubbles are 
interrelated. Thus, the interaction effects between the 
bubbles should be considered in the study of bubble dy- 
namics. Still, rather few modeling efforts have been re- 
ported on the bubble-bubble interaction in cavitation, as 
well as its effects on heat transfer. The present work de- 
scribes a numerical investigation on the dynamics of two 
cavitation bubbles in a heated tube. A three-dimensional 
(3D) model is established to study the interaction effects 
of these bubbles and heat transfer. The liquid and gas 
phases throughout the calculation domain are solved by a 
set of Navier-Stokes equations. The volume of fluid 
(VOF) method is employed to track the interfaces of the 
two bubbles. 


Numerical Method 


In this calculation, the liquid phase is incompressible 
pure water. The bubble is considered to be filled with 
vapor behaving as an ideal gas. The heat transfer between 
the gas and liquid phase is ignored due to a very short 
duration of the bubble evolution. The flow pattern is as- 
sumed to be laminar in the whole domain, which is de- 
scribed by one set of Navier-Stokes equations. 


Governing equations 
The mass conservation equation for the mixture is ex- 
pressed as 


P .¥.(p3)=0 


Ay (1) 
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where p, Vv and ¢ are the fluid density, velocity and 
time, respectively. 
The momentum equation is given by 


<(pi)+¥-(piv) =-vp +v-[u(vi+v¥") + F (2) 


where p is the static pressure, u is the viscosity, F is 
the sum of external body forces and is expressed as 
= PKV a, 


~ l os(p +A) 


where o is the surface tension coefficient, x is the 
curvature of the interface. 
The energy equation is given by the following form 


<(pE)+V-(¥(pE+ p))=V-(kVT,) (4) 


where k is the effective conductivity, Tę is the fluid 
temperature, F is the total energy. 

The liquid-gas interface is calculated by the volume 
fraction equation 


(3) 


0a 


—+v-Va,=0 5 
a 1 (5) 
The volume fraction of the gas phase is given by 
ata,=1 (6) 


where a is the volume fraction, and the subscripts “I” 
and “g” represent the liquid phase and gas phase, respec- 
tively 

The average fluid density and viscosity in each cell are 
calculated according to the volume fraction values of gas 
and liquid, which are expressed as 


P =A + Phe (7) 
H= MQ + MgQy (8) 
The density of the gas phase is given as 
pM 
Pg = (9) 
8 RT; 


where R, is the gas constant, M is the molecular 
weight of the vapor. 

The two interacted bubbles are initially patched in a 
heated tube with the radius of 8 mm, and the calculation 
domain is shown in Fig. 1. The distances from the left 
and right bubble to the centre point of the tube are the 
same (Hiet = Arig). Therefore, behaviours of the two 
bubbles are symmetrical in the calculation. In order to 
eliminate boundary effects, the total length of the tube is 
set to be 160 mm. The grids in the region of -4Rinax<x< 
4Rmax is refined so as to precisely capture the motion of 
the bubbles. It is found that the results are independent 
when the number of grid cells is more than 3.6x10°. A 
constant heat flux of q = 20 W/cm’ is imposed on the 
tube surface. For the pressure outlet boundary, the value 
of Pou = 0.1 MPa is employed. The volume of fluid (VOF) 
method is employed to track the interfaces of two bub- 
bles. The pressure-implicit with splitting of operators 
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(PISO) algorithm is adopted to solve the flow field. The 
pressure staggering option (PRESTO) scheme is used for 
the pressure interpolation at the cell faces. The VOF 
equation is explicitly solved with the geometric recon- 
struction scheme. The second order upwind scheme is 
selected for computing the mass, momentum and energy 
conservation equations. The mathematical model is 
solved with the time step of 1 x 10° s. The convergence 
criterion is 10° for the continuity and momentum equa- 
tions, and 10° for the energy equation. 


L=160 mm 


Men) Height R=8 mm 


Centre point 


Bubbles 


8R max 


Pressure outlet 


yı 
/ 
Z  * Pressure outlet 


Fig. 1 Calculation domain of the present model 
Results and discussion 


The initial radii of the bubbles are Riep.o=Rrignto=0.8 
mm, and values of initial pressure and temperature of the 
bubbles are Piett,0=Pright,o—10 MPa and Thett,0= Tright,o= 1200 
K, respectively. The initial centre distance between two 
bubbles is Hp=3.6 mm, and the initial temperate of the 
liquid in the domain is 7,,=300 K. The numerical results 
are compared with the observation result from the ex- 
periment by Ji et al.'°*!, as shown in Fig. 2. The evolution 
of two cavitation bubbles in a tube obtained by the pre- 
sent simulation agrees well with the experimental pho- 
tography, which proves the accuracy and reliability of the 
model. Due to the interaction effects, dynamic behav- 
iours of each individual bubble are significantly influ- 
enced, which is very different from the case of a single 
bubble near a solid boundary or between the two parallel 
walls. At the beginning, the initially both spherical bub- 
bles grow to a roughly ellipsoidal shape, and the growth 
rate of the inner side of the bubble is lower than that of 
outer side. In the late stages of the growth phase, only a 
very thin liquid film is left between the two bubbles, as 
shown in Fig. 2(a). With further expansion, the two bub- 
bles merge into a bigger one (See Fig. 2(b)). After reach- 
ing the maximum volume, the bubble begins to collapse. 
Figure 2(c) shows that, during the collapse stage, it firstly 
shrinks from the left and right side to the bubble centre. 
Then, the bubble is penetrated by the liquid jets from the 
left and right sides and is separated into two small 
sub-bubbles at the final stage of the first collapse phase. 
As shown in Fig. 2(d), the sub-bubbles are finally pene- 
trated by the liquid jets towards the wall surface and turn 
into a ring shape, which is the second collapse being 
similar to that of a single bubble near a solid boundary or 
between two parallel walls. 
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Fig. 2 Comparison between the numerical and experimental 
results 


Figure 3 shows the typical flow fields around the two 
interacted bubbles at different times. The reference ve- 
locity of the flow jet is represented by the arrow at the 
upper right corner of each drag. It can be seen that at t = 
0.4 ms, in the growth phase of the bubbles, the liquid 
between the bubbles is squeezed out leaving only a very 
thin film. Meanwhile, the outward flow jets with the 
maximum velocity of 4.5 m/s can be observed around the 
bubbles, as shown in Fig. 3(a). After the fusion of two 
bubbles, the merged bubble continually grows until 
reaching the maximum volume. Also, it is shown that the 
speed of the outflow jet gradually decreases during the 
late stage of the growth phase being about only 2 m/s at t 
= 1.4 ms, which is attributed to the weakened effects of 
the bubble expansion (see Fig. 3(b)). During the first 
collapse phase, as can be seen, the opposing inflow jets 
near the left and right sides of the bubble are formed due 
to the contraction of the bubble interface. Meanwhile, the 
liquid micro-jets develop on the upper and lower sides of 
the bubble, which flow outward along the tube after im- 
pinging on the rigid surface, as shown in Fig. 3(c). It is 
shown in Fig. 3(d) that the micro-jets directing towards 
the rigid surface finally penetrate the sub-bubbles in the 
second collapse phase and the maximum velocity of the 
jets increases to about 12 m/s. However, the velocity is 
still much lower than the case of a single bubble near a 
solid boundary or between two parallel walls. This is 
because of the hinder effect of the bubble. In addition, 
the flow fields around the bubbles show that, in the first 
and second collapse phases, stagnation points exist on the 
middle of the tube surface due to the impinging of the 
micro-jets (vortex flow). 

Figure 4 shows the pressure contours around the bub- 
bles with various times. The initial pressure inside the 
bubble is up to 10 MPa, which releases extremely high 
energy at the very beginning of the growth and leading to 
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Fig. 3 Flow fields around the bubbles 


the bubble expansion. However, the motion of bubble 
interface lags behind the transmission of pressure in liq- 
uid. Thus, it can be observed that pressure in the bubble 
is even lower than it in the surrounding liquid, in spite of 
the bubble growth (see Fig. 4(a)~(b)). As time goes on, 
the bubble begins to collapse after growing to the maxi- 
mum volume. When ¢ = 2.8 ms, it can be seen that the 
higher pressure around the bubble lead to the bubble 
contraction, as shown in Fig. 4(c). Thus, inflow jets to- 
wards the tube centre are developed which causes the 
shrinking and splitting of the bubble in the first collapse 
phase. In contrast, the pressure field in the tube centre 
region between two sub-bubbles shows a significant in- 
crease in the second collapse, as shown in Fig. 4(d). 
When ¢ = 4 ms, a maximum value of 1.32 MPa is ob- 
served in the high pressure region, which is the main 
reason for a higher velocity of the impinging jets ob- 
served in Fig. 3(d). 

The deformation of the bubbles and the resulting liq- 
uid flow would affect heat transfer between the liquid 
and the wall. Figure 5 shows the temperature distribu- 
tions of the wall surface along the x direction. On the 
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whole, temperature of the wall surface gradually in- 
creases as the heating time. During the growth phase of 
the bubbles, temperature along the tube surface is almost 
unaffected due to the limited effects of the liquid flow. In 
contrast, the bubble interface gets closer to the wall sur- 
face as well as the resulting liquid flow in the first col- 
lapse phase. Thus, heat transfer in the centre region of the 
wall is enhanced owing to the impinging effects of the 
vortex flow induced by the bubble motion. 

Heat transfer deterioration is found around the centre 
region of the wall, which is mainly attributed to the 
thermal resistance induced by the bubble. During the late 
stages of the second collapse phase, the heat transfer en- 
hanced region is narrowed (in the range of -0.002<x< 
0.002, when 4 ms) due to the limited width of the im- 


pinging jet. 
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Fig. 4 Pressure distributions around the bubbles 
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Fig.5 Distributions of the wall temperature 


Figure 6 displays the changes of the wall temperatures 
with time. The Tw», Two and Ty, max are the wall tempera- 
ture at the far field, the wall temperature at the centre 
point (stagnation point) and the maximum wall tempera- 
ture, respectively. It can be seen that the Tw increases 
lineally under the constant heat flux without any distur- 
bance. The increasing rate of Tw, is much lower than Tw,» 
especially during the late stages of the second collapse 
phase. It is believed that the impinging micro-jet formed 
during the second bubble collapse plays a major role for 
the better effects of heat transfer enhancement. The Ty, 
max 18 almost the same as Tw» before 2 ms. Since then, it 
rises dramatically because the bubble are getting closer to 
the wall surface and therefore induce extra heat transfer 
resistance. For the application of cavitating flow in heat 
transfer, the micro-jet should be strengthened and the 
bubble interface should keep away from the heating sur- 
face by the control of the bubble dynamic behaviours. 
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Fig.6 Variations of the wall temperatures with time 
Conclusions 


On the basis of the VOF method, the interaction of 
two cavitation bubbles in a tube is numerically investi- 
gated together with its effect on heat transfer by using a 
three-dimensional model. The evolution of bubble profile 
results show that the bubble mainly grows along the 
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length (x direction) of the tube because of the wall re- 
striction. The bubble changes to a roughly ellipsoidal and 
only a very thin liquid film is retained between the bub- 
bles. With further expansion, the bubbles finally merged 
into a single one. The outward flow has been observed in 
the bubble growth process. After reaching the maximum 
volume, the single bubble begins to collapse. It firstly 
shrinks from the left and right side to the centre. Then, 
the liquid jets from the left and right sides will penetrate 
the bubble and separate it into two sub-bubbles. The 
sub-bubbles continually shrink and are penetrated by 
liquid jets and turn into a ring shape, which is the second 
collapse being similar with the case of a single bubble 
near a solid boundary or between two parallel walls. Heat 
transfer at the centre region of the wall surface is en- 
hanced owing to the vortex flow and the micro-jet in- 
duced by the bubble contraction and collapse. In contrast, 
the highest surface temperature appears in the surround- 
ing region, which is mainly attributed to the thermal re- 
sistance induced by the bubble. During the late stages of 
the second collapse phase, the heat transfer enhancement 
effect is better due to the impinging of the micro-jet. 
Meanwhile, the value of Tw» has a significant increase 
after 2 ms because of the heat transfer resistance induced 
by the bubble interface. Thus, the bubble interface should 
keep away from the heating surface to avoid the heat 
transfer deterioration. In the future work, the effects of 
such parameters as the initial bubble radius, pressure and 
position on bubble dynamics and heat transfer will be 
discussed in details. 
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